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1.0  Introduction 

Proposed  remedial  acuvities  of  the  shallow  Biscayne  aquiter  at  the  Old  South  Dade  LandtiU  (OSDI.) 
include  extraction  of  ammonia-impacted  groundwater,  treatment,  and  re-injection  ot  the  treated 
water.    Extracuon  of  groundwater  will  occur  from  the  shallow,  G-ll  (freshwater)  aquiter  and  re- 
injccnon  of  treated  groundwater  with  ammoma  concentrations  below  0.5  mg/1  will  occur  in  the 
deeper,  G-lII  (saltwater)  aquifer.   Based  on  the  review  of  the  Remedial  Action  Plan  dated  |ulv,  1998, 
the  National  Park  Scrv^ice  (NPS)  has  expressed  concerns  regarding  the  potential  tor  low 
concenaation  ammonia  discharge  to  Biscavne  Bay  associated  with  migration  of  re  injected  treated 
groundwater.    In  their  letter  dated  January-  25,  1999,  the  NPS  requested:  1)  installadon  of  additional 
groundwater  monitoring  wells  in  the  Bay  ven'  near  to  shore,  and  2)  additional  numerical  mc)deling  to 
further  assess  the  hvdrogeological  interactions  of  the  proposed  remedial  system  near  the  Bay. 

Brown  and  Caldwell  proposed  to  address  the  NPS  concerns  with  additional  Site  characterization  and 
groundwater  flow  modeling  to  support  the  siting  of  a  near-shore  monitonng  location.  Establishing  a 
monitoring  location  to  collect  data  regarding  current  discharge  concentrations  to  the  Bay  will 
provide  a  basis  tor  assessing  potential  discharge  due  to  injection  of  treated  water. 

This  technical  memorandum  summarizes  the  results  ot  additional  Site  characterization  and 
groundwater  flow  modeUng.   The  revised  conceptual  model,  results  of  modeling  activiues,  and 
recommendations  for  siting  a  near-shore  monitoring  location  are  presented. 

A  consensus  meeting  to  discuss  well  siting,  well  mstallation  and  sampling  protocol  will  be  held  prior 
to  initiating  a  water  quality  sampling  event. 


2.0   Conceptual  Model 

The  existing  groundwater  model  (Remedial  AcUon  Plan,  1998)  incorporated  regional  flow  directions 
and  gradients,  canals,  recharge,  evapotranspuauon,  and  groundwater  withdrawals  from  the  opcranng 
landtill  to  the  north  of  the  OSDL  Site  (Figure  1).   To  expand  the  concepmal  model,  two  additional 
components  were  added  to  the  concepruabzation:   discharge  of  groundwater  to  Biscayne  Bay,  the 
saltwater/treshwater  intertace,  and  the  densit\'-dependent  flow  regimes  of  the  near-shore  zone. 

Discharge  to  Biscayne  Bay.   The  issue  of  coastal  freshwater  discharge  to  Biscayne  Bay  is  currently 
the  sub|ect  ot  extensive  research.    Freshwater  discharge  reductions  to  the  Bay  caused  by  the 
lowering  ot  the  water  table  via  managed  canal  and  control  strucnares  has  modified  the  groundwater 
flow  patters  in  southern  Florida.    Quantification  of  groundwater  discharge  to  Bisca\ne  Bay,  not 
directly  measurable  in  the  tleld,  is  part  of  an  on-going  groundwater  modeling  project  bv  the  Lnited 
States  Geological  Sur^-'ey  (LJSGS)  as  part  of  the  South  Florida  Fxosystem  Program. 

Based  on  the  current  understanding,  freshwater  discharge  to  the  Bay  was  obser\ed  and  documented 
prior  to  i960,    breshwater  discharge  was  estimated  to  extend  approximately  MH)  feet  offshore,  based 
on  an  analvsis  ot  freshwater  equipotcntial  (Kohout,  1964).    Since  that  stud\-,  the  general  lowering  of 
tlie  water  table  has  reduced  the  pressure  gradient,  the  ma)or  mechanism  for  freshwater  discharge  to 
the  Bay.    I'reliminarv  data  trom  the  I'SCiS  research  on  discharge  into  the  Bav  support  the  concept  of 
near-shore  discharge;  data  suggest  that  most  discharge  is  limited  to  within  15(1  feet  of  the  shoreline 
(l^angevin,  1999). 
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Saltwater  Interface.    A  result  of  decreased  freshwater  discharge  is  the  intrusion  ot  salr^vater  into 
the  Biscavne  aquifer  system.   A  wedge  of  saltwater  e.xjsts  in  the  aquiter  svstem  undcrK  ing  the 
freshwater  m  the  coastal  region,  and  extends  for  miles  inland.    Data  trom  199S  anahzed  by 
Sosenshein  (1996)  indicates  that  the  toe  of  the  saltwater  zone  is  approximately  two  miles  west  of  the 
Old  South  Dade  Landfill  Site. 

The  contact  between  freshwater  and  saltwater,  termed  the  saltwater  interface,  is  detined  as  the 
10,000  mg/L  Total  Dissolved  Solids  (TDS)  contour  line.   Although  defined  as  a  sharp  contact,  the 
interface  is  acmalJv  a  large  diffuse  zone  of  varying  concentraaons  (Sosenshein,  1996),  formed  bv  the 
mLNing  of  freshwater  with  denser  salr^vater  (Figure  2).  The  mixed  water  zone  is  called  the  zone  ot 
diffusion.     Hydrogeologic  parameters  of  the  aquifer  svstem,  such  as  gradient,  hydraulic  conductivit\- 
and  water  level  flucmaaons,  will  impact  the  thickness  of  this  zone.   Addidonallv,  it  the  aquiter  is 
subject  to  water  level  variations  caused  by  ades,  the  zone  of  mLxed  water  will  be  enlarged  (Fetter, 
1994).   The  resuldng  groundwater  system  is  compnsed  of  two  separate  miscible  fluids:    freshwater 
and  saltwater.   The  densm^  differences  between  the  water  types  exerts  control  on  the  tlow  regime, 
discussed  in  the  following  section. 

Density-Dependent  Flow  Regimes.  In  the  near  shore  aquifer  system  at  Biscayne  Bav,  saLiniacs 
var}^  over  short  distances  both  laterally  and  verdcaUy.  Densit)'  differences  caused  bv  the  varying 
saHnides  separate  the  aquifer  system  into  three  separate,  interconnected  tlow  regimes.    Figure  3  is 
the  locadon  of  the  treshwater/saltwater  interface  based  on  the  most  recent  water  quality  analvses 
from  the  on-site  monitormg  well  network  (OSl  through  OSS)  and  from  saLimt)'  profdcs  completed 
in  the  injecdon  pilot  test  well  (lW-5).    Based  on  the  salinit)'  profiles,  the  zone  of  diffusion  was  a 
mmunum  of  40  feet  thick  (Figure  4).    The  cur\'e  of  the  freshwater/saltwater  interlace  was 
extrapolated  to  the  shoreline,  to  aid  in  the  de\clopment  of  a  conceptual  model  ot  the  saltwater 
interface  and  the  zone  of  ditfusion  in  the  vicimt)-  ot  the  OSDL  Site  and  Biscavne  Bay. 

Figure  5  is  a  west/east  cross  section  of  the  tlc^w  patterns  in  the  freshwater,  diffusion  and  saltwater 
zones.   This  conceptualizatif>n  combines  on-Site  geologic  and  h\  drogeologic  data,  research  on  the 
relationship  between  freshwater  and  salnvater  flow  and  coastal  discharge  in  the  Miami  area 
(Robinson  and  Gallagher,  1999;  Cooper,  1964;  Kohcnit,  1964),  and  preliminar\-  data  from  the  on- 
going L'SGS  modeling  smdv  ot  groundw  atcr  discharge  to  Biscavne  Bav. 

Separate  flow  regimes  are  depicted  tor  freshwater  and  saltwater,  which  are  separated  by  the  zone  of 
diffusion.      The  edges,  or  boundaries,  of  the  zone  of  diffusion  are  restrictions  to  flow.     1  he  upper 
boundary  acts  as  a  barrier  to  flow  of  pure  saltwater,  the  lower  hoiindan'  acts  as  a  barrier  to  flow  of 
pure  freshwater. 

Saltwater  moves  inland  beneath  the  freshwater  lens  in  a  circulatory-  pattern.    Inland  flow  (jf  saltwater 
is  largelv  in  the  deeper  portions  of  the  ac]uifer  (Cooper,  1964),  although  mixing  at  the  interface 
causes  a  component  ot  horizontal  tlow  along  the  length  ot  the  interface. 

Ciroundwater  in  the  freshwater  portion  of  the  aquiter  moves  seaward,  and  encounters  saline  water  at 
the  saltwater/freshwater  interface.    Mixing  occurs  at  the  interface,  which  is  enhanced  In  iitlal 
fluctuanoiis.    Both  the  freshwater  and  mixed  waters  are  less  dense  than  pure  saltwater,  and  mo\e  up 
along  pathlines  indicated  on  ligure  ?>.  Due  to  the  densif\'  gradients,  continuous  mixing  occurs  m  the 
diffuse  zone  as  water  mo\  es  upward  along  the  interface.     lc)\vard  the  upper  j^jortion  of  the  aeiuifer. 
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waters  m  the  zone  of  diffusion  are  representanve  of  combined  treshwater  and  saltwater  tioni  all 
laxcrs  in  the  aquifer. 

Discharge  into  the  Bay  consists  of  freshwater  from  the  shallow  aquifer  system  and  mixed  water 
from  the  zone  of  diffusion  (see  enlargement.  Figure  5).    Most  groundwater  discharge  was  estimated 
to  occur  within  150  feet  of  shoreline.  Based  on  research  in  the  vicinih'  (Langevin,  1999),  a 
freshwater  component  extends  offshore  approximately  91)0  feet.   This  represents  a  limit  tor  the 
lateral  extent  of  discharge,  and  it  is  assumed  that  discharge  with  any  freshwater  component  at  all 
occurs  witliin  this  near-shore  zone. 

Impacts  of  Injection.   As  proposed  in  the  R^\P,  treated  groundwater  from  the  shallow,  G-Il 
aquifer  (freshwater)  will  be  re-mjected  into  the  deeper,  G-III  aquiter  (saltwater)  for  disposal.    The 
injection  intenal  is  from  -60  to  -100   ft  below  sea  level  (bsl).  The  injecuon  wells  are  located 
approximately  4,500  feet  from  the  shoreline  (see  in)ecuon  well  IW-5  on  Figure  6).  The  in)ecdon 
interval  at  this  locadon  is  within  the  diffusion  and  saltwater  zones.  Freshwater  mjected  into  these 
zones  will  produce  a  freshwater  mound  that  rruxes  and  moves  up  until  reacliing  a  densit}' 
equilibrium.  As  the  freshwater  moves  up,  the  densirv'  gradient  becomes  smaller,  reducing  the 
mechamsm  for  vertical  movement.   As  the  driving  mechamsm  is  diminished,  horizontal  gradients 
will  begin  to  influence  flow,  and  a  horizontal  component  of  flow  will  cause  injcctate  to  move  up  and 
east  toward  the  Bay.    Injectate  will  likely  reach  a  densir\'  equilibrium  near  the  upper  boundar\'  of  the 
zone  of  diffusion,  just  below  the  freshwater  zone  (Figure  6).  Continuous  rruxing  with  saltwater  and 
freshwater  will  occur  as  the  injectate  moves  along  the  boundary.   The  semi-con tunng  layer  is  a 
barrier  to  upward  migradon  of  water  mto  the  shallow  zone,  as  indicated  by  the  small  verdcal  arrows 
in  that  region. 


3.0   Numerical  Groundwater  Model 

The  expanded  conceptvial  model  incoiporates  some  complexities  that  were  not  contained  in  the 
active  modeled  area  of  the  existing  groundwater  flow  model  tor  the  Site.   These  include  the 
freshwater/saltwater  interface,  separate  tlcnv  regimes  for  treshwater  and  saltwater,  treshwater 
discharge  to  the  Bay,  and  groundwater  flow  in  the  aquifer  underlying  the  Bay.    A  review  of  modeling 
studies  ot  analrjgous  systems  was  conducted  to  assess  the  difterent  simulation  methods  and  their 
relative  successes.    To  best  simulate  the  densir\-dependent  interactions  between  the  treshwater  and 
saltwater  systems,  and  to  achie\'e  the  resolution  required  to  depict  shoreline  discharge,  tlnite  element 
models  utilizing  pressure  head  and  sahnir^'  in  the  governing  equaUons  for  groundwater  flow  and 
contaminant  transport  were  most  prevalent. 

MODFLOW  (McDonald  and  Flarbaugh,  1988),  a  three-dimensional,  flnite  difference  groundwater 
model  that  calculates  groundwater  heads  and  budgets,  was  used  to  develop  the  existing  model. 
MODIT.OW  IS  ncjt  densm-dependent,  and  its  finite  difference  structure  is  less  applicable  to  non- 
linear boundary  conditions.   This  precluded  direct  modeling  of  the  dual-densirv  system  with  the 
existing  model.  Model  development  therefore  focused  on  the  freshwater  flow  regime. 

the  actne  area  ot  the  existing  numerical  groundwater  model  was  expandcLl  to  include  the  shoreline 
and  near-shore  regions  in  Bisca\ne  Bay.     The  existing  fmire  difference  grid  ranged  frf)m  60  foot 
squares  on  sire  to  2l)(l-t()(>t  squares  toward  the  edges  of  the  modeled  area.    Tlie  refined  area  of 
smaller  6()-toot  grid  cells  was  expanded  to  include  the  shoreline  and  near-shore  regions.    Sources  of 
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water  into  the  expanded  groundwater  system  include  leakage  from  canals,  intiltration  ot 
precipitation,  and  underflow.   Water  leaves  the  groundwater  system  via  discharge  into  Biscavne  Hav, 
leakage  into  canals,  evapotranspiration,  and  underflow. 

The  observed  hydrogeologic  components  of  the  remedial  system  and  the  aquiter  system  are 
discussed  below,  followed  by  the  approach  taken  to  simulate  each  component  in  the  numerical 
model. 

The  Remedial  System.    The  remedial  system  at  the  OSDL  will  consist  of  extraction  wells,  water 
treatment,  and  in)ecrion  wells  to  dispose  of  treated  groundwater.  Five  extraction  wells  pimipmg  at 
rates  that  vary  from  101  to  156  gallons  per  minute  (gpm)  will  hydraulically  control  the  shallow 
groundwater  system,  capmrmg  ammonia  rmpacted  groundwater  underlying  Area  III  of  the  OSDL 
(Figure  7).   Total  dailv  treatment  volumes  wlU  be  approximately  925,000  gallons.   Once  treated, 
groundwater  wUl  be  re-mjected  into  the  deeper  (G-III)  aquifer  of  the  Fort  Thompson  Formation, 
below  the  10,000  TDS  concentration  contour  delineating  the  boundary  between  freshwater  and 
saln.vater.  A  total  of  five  injection  wells  each  injecung  127  gpm  to  dispose  of  the  treated 
groundwater  wiU  be  located  approximately  300  feet  east  (downgradient)  of  the  pumping  wells 
(Figure  7).    The  m)ection  wells  will  be  surrounded  bv  percolation  ponds  which  dispose  ot 
stormwater  collected  as  runoff  from  Area  111.    During  large-scale  storm  events,  runoft  will  drain  into 
the  injecdon  wells. 

Injecdon  of  treated,  fresh  groundwater  below  the  10,000  TDS  concentradon  contour  will  change  the 
vertical  gradients  in  the  vicrmn^  of  the  injection  wells.   An  upward  vcrdcal  gradient  currendy  exists, 
which  will  be  enhanced  by  1)  the  posidve  pressure  response  in  the  deeper  system  due  to  injection, 
and  2)  the  introduction  of  freshwater  into  the  saline  envu-onment,  creating  dcnsitv  diflercnccs  that 
increase  upward  movement.    Mounding  due  to  mjection  was  simulated  (Remedial  Action  Plan, 
1998),  was  approximately  0.17  feet  (0.17  feet),  and  was  restricted  to  the  immediate  vicinit)'  of  the 
well.  UpweUing  due  to  density  differences  cannot  be  directlv  simulated  with  the  MOD  FLOW  model. 

Pumping  and  extracdon  were  simulated  with  wells  in  the  groundwater  model.    Percolarion  ponds 
were  sunulated  using  a  real  recharge  rates  set  at  the  design  percoladon  rate.    Bv  var-sing  the  well 
injectKm  rates  and  the  duration  of  percolation,  both  daily  operating  concUdons,  small  scale  storm 
events,  and  a  lOO-year  storm  event  were  simulated. 

Biscayne  Aquifer.    The  thickness  of  the  Biscavne  aquifer  in  the  viciniU'  of  the  OSDL  is 
approximately  100  feet,  based  on  subsurface  investigations  conducted  on  Site  (Brown  and  C^aldwell, 
1999).    The  Biscayne  aquifer  in  the  vicimt)'  of  the  site  is  comprised  of  the  Miami  Limestone  and  the 
Fort  '1  hompson  hormauon,  which  are  separated  bv  a  semi-confining  layer  that  occurs  at  the  contact, 
approximately  25  feet  below  sea  level  (ft  bsl).  The  semi-confining  layer  does  not  eliminate  vertical 
flow,  but  acts  as  a  barrier  to  \erncal  groundwater  mo\  emenr. 

Separate  layers  within  the  numerical  model  were  used  to  represent  the  Miami  Limestone,  the  Fort 
Thompson  lormation  and  the  semi-con  fining  laver  at  the  contact.    Flvdrogeologic  parameters  were 
estimated  ba^cd  on  pumping  tests  in  the  shallow  Miami  l.unestone  and  mjecnons  tests  in  the  deep 
fort  ihonipson  formation,  designed  to  assess  the  h\clraulic  connection  between  ac|uifers. 
Relatively  high  values  of  hydraubc  conducuvm'  were  utili/cd  for  the  Miami  Limestone  and  fort 
Thompson  formation  (.\0(KI  to  3,500  ft/dav);  relativeK  low  values  of  hvdraulic  conductivirv  were 
u>ed  tor  the  semi-conflning  layer  (50  ft/day).    X'erfical  leakage  in  the  semi-conflning  laver  was  set 
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nvo  orders  of  magnitude  lower  than  in  the  Miami  Limestone  and  Fort  Thompson  i-ormation.    Due 
to  the  large  disparm-  in  conductiviues,  preferential  flowpaths  are  pnmarilv  horizontal  and  limit 
hydraulic  communicauon  between  aquitcr  units. 

Flow  Regime.  Groundwater  flow  across  the  OSDL  is  to  the  east/southeast.  The  flow  regime  in  the 
vicinih^  of  the  Site  is  governed  by  two  major  hvdrogeologic  controls:   diurnal  tides  and  canals.    The 
Goulds  canal  runs  west/east  along  the  northern  boundan'  of  the  entire  Site,  a  spillover  gate  drains 
water  from  the  canal  into  the  Bay.   The  L31E  canal  connects  to  the  Cioulds  Canal  and  runs 
north/south,  between  Areas  11  and  111  (Figure  1).   The  L31E  is  an  uncontrolled  flood  management 
strucmre.   The  canals  are  excavated  into  samrated  Miami  Limestone  and  are  assumed  to  have  direct 
hydraulic  connection  to  the  aquifer.   Two  large  west/east  canals  with  control  structures  are  located 
to  the  north  (Black  Creek  canal)  and  south  (Princeton  canal)  of  the  Site.   Canal  levels  are  generally 
lower  m  the  summer  m  anticipation  of  large  influxes  of  water,  and  higher  in  the  winter. 

Canals  were  simulated  with  the  River  Package  of  MODFLOW,  an  internal  boundan'  that  allows 
water  to  move  mto,  or  out  of,  the  groundwater  system,  depending  upon  the  gradient  between  the 
nver  stage  and  the  aquifer. 

Two  disnnct  seasons  are  recogruzed  m  this  coastal  climate:   the  dry'  season  representing  low-water 
months  during  the  winter  and  spring,  and  the  wet  season  representmg  summer  and  fall.   Sevenr\'-five 
percent  of  the  total  estimated  rainfall  of  55  inches  per  year  falls  during  the  summer  months  (Merritt, 
1996).    In  addition  to  seasonal  water  level  variations,  the  aquifer  system  m  the  vicinirv  of  the  site  is 
subject  to  tidal  fluctuations.   Momtonng  of  water  level  elevations  during  a  12-hour  tidal  cycle 
showed  that  adal  vanadons  were  dampened  by  canal  unpacts.   Areas  I  and  II  (Figure  7)  showed  a 
ndal  response  of  approximately  0.5  ft  and  Areas  111  and  IV,  west  of  the  L31E  canal  showed  a  tidal 
response  of  approximately  0.1  ft.    Flow  durecticju  and  gradient  were  mamtained  throughout  the 
diurnal  perturbadons. 

The  numerical  model  was  designed  to  simulate  the  nominal  system,  midway  bet^veen  the  high  and 
low  water  season.  I'ldal  vanations  were  not  djiectlv  simulated,  water  levels  reflected  average  tidal 
conditions. 

Saltwater  Interface.   The  lower  boundar)-  of  the  zone  of  diffusion  was  conceptualized  as  a  barrier 
to  freshwater  flow  (Figure  5)  that  forms  a  sharp  interface  between  the  freshwater  and  salr^vater 
regimes.   Although  the  zone  of  diffusion  is  know  to  be  a  transient  feature,  not  well  represented  by  a 
shaip  bf)undar\',  the  simplifying  assumption  of  a  shaip  boundar\'  was  made  for  simulation  puiposes. 
fhis  representation  of  a  sh;iip  boundary  is  consenatue  because  the  discharge  is  restricted  to  flow- 
along  the  interface  and  not  dispersed  in  a  rriLxing  zone  rhrf)ughout  the  interface. 

I  he  boundary'  separating  treshwater  from  saltwater  was  simulated  with  horizontal  flow  barriers  in 
MODl'LOW.    Horizontal  flow  barriers  are  located  between  cells  and  use  a  barner  transmissivity  to 
lower  the  hori;^ontal  branch  conductance  between  the  rvvo  cells  that  the\'  separate.    Barriers  to  flow 
were  positioned  vertically  along  the  saltwater  front,  and  staggered  between  layers  to  approximate  the 
cun mg  boundan-  condition. 

Sources/Sinks.    Water  eiiters  the  acti\-e  modeled  :irea  as  underflow  from  upgradieiit  sources,  canal 
leakage,  and  intilrrarK>n  froiii  precipitation.    Water  leaves  the  actne  mr. deled  area  :is  underflov 
discharge  into  the  Bay,  leakage  into  the  canals,  and  e\  apotranspirati 
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L'nderflow  in  and  out  of  the  active  model  area  was  simulated  with  constant  head  boundanes  to  the 
west  and  east.   The  western  constant  head  boundar)'  was  set  at  2  teet,  and  the  eastern  boundar\', 
moved  approximately  one  mile  offshore  during  model  expansion,  was  set  at  sea  level.   Canal  leakage, 
discussed  previously,  was  simulated  with  River  cells.    Impacts  ot  precipitation  and 
evapotranspiration  were  simulated  with  the  Recharge  and  Hvapotranspiration  Packages  of 
MODFLOW,  respectively.   Groundwater  discharged  mto  the  Bay  leaves  the  groundwater  s\stern, 
and  moves  into  the  surface  water  system.  Conceptually,  this  discharge  exits  the  groundwater  model 
domain.  Shoreline  discharge  was  simulated  with  the  Dram  Package.    Drains  are  interior  model 
bc:>undaries  that  remove  water  from  the  modeled  area  when  water  levels  are  higher  than  a  user- 
defined  stage  in  the  dram  cell.  PreUminan'  results  from  the  L'SGS  shoreline  discharge  mt)del  were 
not  yet  available  for  calibration  to  actual  flux  rates;  rough  estimates  of  discharge  were  calculated 
witlim  the  model,  based  on  a  drain  stage  ot  ().5ft  msl.     Figure  8  shows  all  boundar\'  conditions  used 
for  the  expanded  groundwater  model. 


4.0  Results 

L'se  of  the  honzontal  flow  boundan'  successfully  simulated  the  vertical  components  of  freshwater 
flow  depicted  in  the  concepmal  model.    Model  results  incoiporated  a  component  of  discharge  into 
the  Bay  designed  to  remove  water  from  the  groundwater  system  near  the  shoreUne.  Freshwater 
discharge  to  the  Bay  ranged  from  60  to  80  ft/day  per  foot  of  shorelme.  Due  to  msufficient  data 
regarding  the  obser\'ed  geometn-  of  the  saltwater  mterface  at  the  shoreline,  the  simulated  discharge 
was  not  calibrated  to  field  conditions.   Although  discharge  was  an  uncalibrated  component  of  the 
acjuifer  system,  based  on  background  research  and  our  understanding  of  the  acjuifer  system, 
discharge  rates  in  this  range  are  plausible. 

Particle  tracking  simulations  to  calculate  flow-paths  of  groundwater  were  run  with  the  expanded 
model.    Particle  tracking  simulations  are  based  on  the  ceU-by-cell  fluxes  and  directions  of  flow  that 
are  calculated  m  the  MODFLOW  groundwater  flow  model.    These  outputs  are  used  by 
MODPA'FFI  (IV^Uock,  1994)  to  calculate  velocit)'  and  direction  of  flow  of  a  particle  "released"  into 
the  flow  system.    The  MODPA  I'll  module  was  used  to  track  groundwater  injected  between  60  and 
100  feet  below  sea  le\el  (ft  bsl).  Figure  9  shows  the  flow-paths  of  water  particles  released  in  the 
mjection  interv^al.  The  arrows  on  the  flow-path  lines  are  one-year  inter\'als.   Released  particles 
reached  the  shoreline  after  approximately  4  to  6  years.  Note  that  the  estmiated  time  to  move  up 
through  the  aquifer  system  and  reach  the  shorehne  is  based  on  the  nominal  flow  system,  midway 
between  high  and  low  tides.    Tidal  impacts  and  ma)or  fluctuations  in  water  levels  associated  with 
seasonal  recharge  are  not  reflected  m  this  estimate.     Further,  because  the  interface  is  represented  as 
a  sharp  mterface,  the  simulated  flow-paths  are  likely  conservative  and  represent  the  lower  bound  of 
time  to  discharge  into  the  Ba\ . 


5.0  Recommendations 

Well  Site,    .\ddirional  background  iinesrigations  were  conducted  to  idenrif\-  related  stuilies  and 
existing  near-shore  monironng  locations  to  maximi/e  the  benchts  from  the  proposed  new  ucll.    (  )n- 
gf>ing  inyestigarions  ixlared  to  groundwater  or  surface  water  are  l)eing  conducted  1)\  the  I  'uited 
States  Geological  Suney  (L'SGS),  the  National  (  )ceanic  and  .\fmospheric  Admiiustration  (N(  )A.\), 
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the  National  Park  Scn-ice  (NPS),  the  Armv  Coq^s  of  Engineers  (COE),  and  l-lorida  Inrcrnarional 
Universit}-  (FIU).  Contacts  were  made  with  these  agencies  to  locate  existuig  offshore  monitoring 
wells  in  the  vicimt}-  of  the  Site. 

The  USGS  well  cluster  adjacent  to  Black  Pomt  (15P1A  and  BPIB)  is  closest  to  the  OSDL  Site.  The 
wells  are  located  approximately  2,U()()  feet  offshore,  too  far  to  be  useful  in  the  assessment  ot 
potential  concentrations  of  ammonia  in  near-shore  Bay  discharge,  but  located  at  what  could  be  tlie 
endpoint  of  a  groundwater  monitoring  transect  of  wells  across  the  Site  and  into  the  Bay,  once  the 
near-shore  monitoring  well  is  installed. 

Based  on  the  conceptual  and  numerical  groundwater  models,  the  optimal  location  tor  the  proposed 
new  near-shore  monitoring  well  is  adjacent  to  the  OSDL  Site,  witlnn  50  feet  ot  the  shoreline,  as 
shown  on  Figure  10.   This  location  is: 

■  midwav  ben-veen  the  Goulds  Canal  to  the  north  and  the  Princeton  canal  to  the  south,  reducing 
potential  influences  of  canal  discharge  on  water  qualitv'  sampling  events; 

■  inland  of  the  USGS  Black  Point  monitoring  wells  in  the  Bay,  forming  a  transect  of  monitoring 
wells  across  the  Site  and  out  into  the  Bay;  and 

■  immediately  downgradient  of  the  landfill  and  within  the  groundwater  discharge  zone,  thus 
representative  of  the  concentrations  of  ammonia  that  will  discharge  into  the  Bay. 

However,  the  well  location  is  within  the  mangrove  frmge  present  along  the  shoreline  of  the  Bay. 
The  mangrove  system  has  been  identitled  as  a  namral  producer  of  ammonia  (NPS,  1999),  and  the 
water  quality'  results  from  this  location  may  be  impacted.    It  is  requested  that  the  NPS  aid  in  the 
quantitication  of  the  potential  contribution  of  the  mangrove  tringe  to  the  water  qualit\'  results  to 
better  assess  nitrogen  speciauon  in  the  system. 

Sampling  Protocol.   To  coincide  with  existing  on-Site  monitoring  and  take  advantage  of  available 
data  from  the  nearby  monitoring  wells,  the  new  proposed  well  should  be  a  multi-port  well,  or  a  well 
cluster,  to  provide  a  verdcal  protile  of  nitrogen  species.      Suggested  sampling  intervals  are:   surface 
water,  1-6  ft  below  the  rock  surface  (to  be  consistent  with  the  existing  off-shore  wells);  15-20  ft 
below  the  rock  surface  (to  be  consistent  with  the  existing  off-shore  wells);  and  40  to  50  ft  below  the 
rock  surface,  to  assess  cc^ncentradons  in  the  deeper  G-lll  aquifer. 

A  sampling  event  to  include  all  on-Site  and  off-shore  monitoring  wells  is  recommended,  to  provide 
a  "snap  shot"  of  concentrations,  assess  the  dominant  geochemical  processes  at  work  in  the  aquifer 
systems  and  improve  our  understanding  of  the  subsurface  nitrogen  speciation  downgradient  of  the 
landtill.  To  accomplish  this,  samples  must  be  anaUv.ed  for  the  full  suite  of  nitrogen  species.    It) 
more  closely  tie  in  with  existing  monitoring  projects,  it  is  recommended  that  split  samples  be 
collected  and  analyzed  by  the  respective  labs  for  Miami-Dade  Solid  Waste,  the  CSGS,  and  ML'. 
These  three  entities  have  existing  databases  of  water  qualm-  analyses  and  slightly  different 
methodologies  for  sample  collection  and  presenation.     Split  samples  will  remove  some  of  the 
ambiguities  inherent  in  comparing  results  using  differing  methodologies.    .\  consensus  meeting  to 
discuss  the  well  installation  and  define  a  sampling  plan  will  be  held  subsequent  to  the  submittal  of 
this  memorandum. 
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6.0  Summary 

■  ConccprualJy,  the  aquifer  svstem  near  tlie  Bay  is  density-  dependent.   SimpUhm^  assumptions 
were  made  to  simulate  key  components  ot  the  freshwater  system. 

■  Expanded  conceptual  and  numerical  models  of  the  OSDL  Site  were  utilized  to  assess   the 
hydrogeologic  components  of  the  shoreline  environment.  Revisions  to  the  groundwater  model 
were  successful  m  simulating  the  vertical  components  of  treshwater  flow  and  shoreline  discharge 
depicted  in  the  conceptual  model.    L'ncaLibrated  shoreline  discharge  ranged  trom  60  to  (SO 

ft  /dav  per  toot  ot  shoreUne.     Particle  tracking  was  performed  to  determine  the  optimal  loauion 
for  the  propf)sed  new  monitormg  well. 

■  Recommendations  on  the  well  design  and  sampling  protocol  were  made,  and  were  intended  to 
take  advantage  of  existing  databases,  minirmze  ambiguides  m  results,  and  provide  an  overview  of 
nitrogen  speciation  along  a  transect  from  the  OSDL  Site  to  die  Bav.   The  recommended  well  site 
IS  withm  50  feet  of  the  shoreline. 

■  The  NFS  is  requested  to  aid  in  the  quandfication  of  the  namral  production  of  ammonia  bv  the 
mangrove  tringe,  in  order  to  more  accurately  assess  the  nitrogen  speciation  across  the  proposed 
monitoring  transect. 
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